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^ i ABSTRACT 

The HERA diffractive structure function data are interpreted in terms of 'diffractive parton 
distributions' which satisfy DGLAP evolution. These distributions are modeled assuming 
that the scattering takes place off a colour singlet 'pomeron' target. A quantitative test of 



X 



the universality of diffractive parton distributions is proposed. 



* Based on talks presented at the Workshop on HERA Physics, Durham, September 1995 
and the Workshop on Deep Inelastic Scattering and Related Phenomena (DIS96), Rome, 
April 1996. 



Measurements of deep inelastic scattering events with a large rapidity gap at HERA [|T], 
0, |3|, |], [|] have generated renewed interest in the idea of 'diffractive hard scattering'. It 
now appears that in a significant fraction of these events there is a diffractively scattered 
proton in the final state [EJ. One can therefore introduce the idea of 'diffractive parton 
distributions', i.e. parton distributions in the hadron under the condition that the hadron is 
diffractively scattered. Just as the total deep inelastic structure function, measured in the 
process 7*p — > X, can be written as a sum over parton distributions times a short-distance 
part, F 2 (x B j,Q 2 ) = Eg fdx f q/p (x,fx 2 ) F 2q (x Bi /x,Q 2 , fi 2 ) with F 2q = e\ 5(l-x B j/x)+0(a s ), 
we may define a similar decomposition for the diffractive structure function, measured in the 
process 7*p — > pX, 

dF 2 D (x B j,Q 2 ;x v ,t) _ ^ r df q/p (x, fi 2 ] x v , t) ^ /aty q2 ^ 2 \ ^ 



dx-pdt q J dx-pdt \ x 

Here 1 — xp (xp <C 1) and t (\t\ <C Q 2 ) are respectively the longitudinal energy frac- 
tion and the t-channel momentum transfer of the proton in the final state. The quantities 
df q/ p (x, /i 2 ; xp,t)/ dx-pdt can then be regarded as diffractive parton distributions.^ A rather 
detailed theoretical discussion of such diffractive structure functions has recently been pre- 
sented in Ref. ||. In particular it has been shown that an operator definition exists, and 
that the diffractive distributions should satisfy the same (DGLAP) evolution equations as 
the usual parton distributions. However it is not yet clear whether the concept of short- 
distance factorization generalizes to any diffractive hard scattering process, for example the 
production of large-mass Drell-Yan lepton pairs or large jets in hadron-hadron colli- 
sions ||. Even if factorization is violated for diffractive hard scattering in such collisions, 
the effect may be weak at high energy. It therefore seems to us not unreasonable to assume 
the approximate validity of universal factorization and to test its consequences. The purpose 
of the present study is to explore further the idea of hard diffractive factorization by using 
the HERA F 2 data to model diffractive parton distributions, and then to use these to make 
quantitative predictions for other diffractive hard-scattering processes, in particular for the 
production of W bosons in pp collisions at the Tevatron. 

An important property of the HERA diffractive events [IJ is the approximate factorization 
of the structure function F 2 D (integrated over t) into a function of xp times a function of 
P — x B - s /xp\ F 2 ~ Xp n J r (P,Q 2 ). This, together with the observed rapidity-gap topology of 
the events, strongly suggests that the deep inelastic scattering takes place off a slow-moving 
colourless target V 'emitted' by the proton, p — ► Vp, and with a fraction xp < 1 of its 
momentum.^ If this emission is described by a universal flux function fp(xp,t)dxpdt, then 
the diffractive structure function F 2 D is simply a product of this and the structure function 
of the colourless object, F^^.Q 2 ). Since the scattering evidently takes place off point- 
like charged objects, we may write the latter as a sum over quark-parton distributions, i.e. 
F 2 (P,Q 2 ) = PJ2 q e l a v(P,Q 2 ), in leading order. In this way we obtain a model for the 

1 Similar quantities ('fracture functions') were introduced in Ref. 0. 

2 This physical picture refers to the infinite momentum frame of the proton. Although it is rather different, 
it is not a priori inconsistent with the so-called aligned quark picture fit! which is valid in the rest frame of 
the proton. 
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diffractive parton distributions introduced above: 



Many recent studies have analysed the HERA diffractive structure function data using 
this theoretical framework. A popular choice is to assume that the colour-neutral target 
is the Regge pomeron, in which case the emission factor f-p is already known from soft 
hadronic physics (for a review see Ref. I n its simplest version, this model would predict 

F 2 D (x,Q 2 ;xp,t) = fp(xp,t) F£((3,Q 2 ) with fp(xp,t) = Fp(t)x 2 p v{ ^~ 1 , i.e. a factorized 
structure function with n pd 2op(0) — 1 ~ 1.16. This model is based on the notion of 'parton 



constituents in the pomeron' first proposed by Ingelman and Schlein [12] and supported by 



data from UA8 fl3| . In such a model, a modest amount of factorization breaking, such as 
that observed in the more recent HI and ZEUS analyses j|, |5|, [|, could be accommodated 
by invoking a sum over Regge trajectories, each with a different intercept and structure 
function: 

F 2 D (f3, Q 2 ; x v , f) = £ F n {t)x 2 ^ {t) - 1 F? (/3, Q 2 ) , (3) 

K 

which would yield an effective n which depends on f3 but is approximately independent of 
Q 2 .f\ Note that since in practice the variables xp and t are integrated over, what is measured 
is a linear combination of F^ structure functions or, equivalently, the parton distributions 
in an effective colour- neutral target: 

/ dx P dtF 2 D (P,Q 2 ;x P ,t)=J2 A n F 2 R (P,Q 2 )=(3j2 e l E ^W(/3,Q 2 ), (4) 

where the coefficients An are independent of (3 and Q 2 . Since the DGLAP equations are 
linear in the parton distributions, the Q 2 evolution of the integrated structure function F 2 D 
should also be calculable perturbativelyQ In the present study we assume, for simplicity, 
pomeron exchange only (i.e. 1Z = V) and use the parametrization of Ref. |L5] for a(t) and 



F-p(t). The x-p dependence of the diffractive structure function predicted by this type of 'soft 
pomeron' model is roughly consistent within errors with the HI ||] and ZEUS 0, |6j data, 
although there is some indication from the latter that a somewhat steeper x-p dependence is 
preferred. 

Various models for the parton distributions in the 'pomeron' f q /p(x,fi 2 ) have been pro- 
posed, ranging from the two extremes of mainly gluons to mainly quarks. Recent studies in 



the framework of perturbative QCD DGLAP evolution can be found in Refs. EL 16, 17, 1 



El, E0|. As we shall demonstrate below, models of both types can be constructed to agree 



with the HERA data. A key issue concerns the existence of a momentum sum rule for the 
colourless object V. This is a matter of some dispute, and there is indeed no theoretical proof 



3 For a recent quantitative study see Ref. |Q. 

4 In principle, for the integrated structure function there is a correction to the evolution equations from the 
small but finite probability that the final-state proton is produced in the hard-scattering process. However if 
the integration is only over a limited region in t, as in the present context, this correction will be completely 
negligible. 
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that such a sum rule should exist. Of course because it is the product of f q /-p and f-p that 
appears in the expression for the structure function, one can simply impose a momentum 
sum rule on the parton distributions and absorb an overall normalization A/", unchanged by 
Q 2 evolution, into fp. This is the approach we shall adopt here. 

Since our purpose is to explore the consequences of a parton interpretation of deep 
inelastic diffractive scattering for other processes, we consider three qualitatively different 
models of diffractive parton distributions. The first is purely quark-like at the chosen starting 
scale /Xq, with gluons generated only at higher scales through standard DGLAP evolution. 
The second contains a mixture of quarks and gluons at /x 2 ,, and the third is predominantly 
gluonic. Since in each case the quark content is constrained by structure function data, and 
since we choose to impose a momentum sum rule, the overall normalization factors M are 
different in the three models (see below). In each case the starting distributions are chosen 
to give satisfactory agreement with the HI and ZEUS || data. Since the errors on these 
data are quite large, it is not necessary to perform detailed multiparameter fits. The starting 
distributions (at /Xq = 2 GeV 2 ) are given in the following Table. 



Model 






N 


1 


0.314x 1 / 3 (l -x) 1 / 3 





1.62 


2 


0.2x(l-x) 


4.8x(l - x) 


2.85 


3 


0.081x(l -x) 05 


9.66x 8 (l -x) 02 


1.57 



Here q refers to an individual light-quark distribution with SU(3) flavour symmetry assumed, 
i.e. q = u = u = d = d = s = s, so that the momentum sum rule constraint at /Xq is 
Jq 1 dx x(6q + g) = 1. Charm quarks are generated by massless DGLAP evolution (g — > cc) 
at higher scales. The Table also lists the normalization factors required to fit the data. We 
emphasize that these are correlated to our use of the Donnachie-Landshoff parametrization 
of f-p(xp,t)^\ Different emission factors will lead to different normalizations. Fig. [I] shows 
the quality of the fits to the HI diffractive structure function data 0. The fits to the ZEUS 
data are comparable in quality and are not shown. More details about the fits and the 
different Q 2 evolution in the three models will be given elsewhere pi . 



The concept of 'universal pomeron structure' can be tested in hard diffractive processes 
in hadron-hadron collisions. The cleanest process to study would appear to be weak boson 
production at the Tevatron pp collider. This was first studied in Ref. [22], where it was 
estimated that the single diffractive component of the total W cross section could be as 



large as 20%. More recently, the CDF collaboration [23| has searched for such diffractive 
events, and derived a preliminary upper limit on the single diffractive cross section of 'a 
few per cent'. Using the three models discussed above, it is straightforward to compute the 
single diffractive W cross section: 



a (W) ~ 2q p ®Aff P q v , 



(5) 



where the factor of two corresponds to either the proton or antiproton being quasi-elastically 
scattered. To avoid uncertainties from high-order corrections it is sensible to normalize this 

We integrate the x-p and t variables over the range appropriate for the HERA experiments. 
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Figure 1: Fits (solid line: Model 1, dashed line: Model 2, long-dashed line: Model 3) to the 



HI Fj> data @. 



prediction to the total W cross section, a (W) ~ q p <E> ? p . For the parton distributions in 
the proton we use the MRS(A') set f24j| . The factorization scale is chosen to be Q = Myy, 
and we sum over (four) flavours of quarks in the initial-state proton and pomeron. Following 
Ref. [22|, we define 'single diffractive' events by x-p < 0.1, integrating over all t. In practice, 
the events are defined by rapidity gaps of a certain minimum size, and therefore the observed 
diffractive cross section must be corrected to the theoretical prediction based on, say, x-p < 



0.1 using a Monte Carlo simulation [23|. With the above choice of parameters and cuts we 
find 

Jaw [ 5-3% Model 1 

6.5% Model 2 (6) 
7.4% Model 3 



a' 



a 



(W) 



The important point to note here is that the model predictions are quite similar: the hy- 
pothesis of hard diffractive factorization has yielded a well-constrained prediction for the 
single diffractive cross section. This result can be understood in terms of the evolution of 
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the quark densities in the pomeron. The single diffractive W cross section at the Tevatron 
samples the quarks in the pomeron at (x q /-p) ~ 0.4. At low Q 2 the quark distributions at 
this x value are constrained by the HERA data to be roughly the same. As Q 2 increases 
the distributions diverge, reflecting the quantitatively different gluon contributions to the 
DGLAP evolution. However at Q 2 ~ 10 4 GeV 2 , the relevant value for W production, the 
difference between the quarks in the three models is still not very large, and the predictions 
for o SD (W) are correspondingly similar, see Fig. l.f] We stress that failure to observe a 
diffractive W cross section of the order of the values given in Eq. (^) would cast serious 
doubt on the 'universal pomeron structure' hypothesis. 

In summary, we have shown how the HERA diffractive structure function data can be 
understood in terms of diffractive parton distributions, which satisfy DGLAP evolution and 
can be modelled in terms of various combinations of quarks and gluons in an effective colour- 
neutral target. We have discussed the concept of the universality of such distributions, and 
shown how the measurement of the single diffractive W cross section at the Tevatron will 
provide a stringent test of the universality property. 
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